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PACS numbers: 71.10.Fd, 71.10.Hf, 71.10.Pm, 71.30.+h 

In their Comment Q], Sandvik, Sengupta, and Camp- 
bell present some numerical evidences to support the ex- 
istence of an extended bond-order- wave (BOW) phase at 
couplings (U, V) weaker than a tricritical point (Ut,Vt) 
0, El in the ground state phase diagram of the one- 
dimensional half-filled U-V Hubbard model. They claim 
that their results do not agree with the phase diagram 
proposed in my Letter Q , which shows a BOW phase for 
couplings stronger than the critical point only. However, 
I argue here that their results are not conclusive and do 
not refute the phase diagram described in the Letter. 

First, while the parameter U — At used in the Com- 
ment is smaller than the tricritical coupling XJ t found in 
Ref. 0, it is larger than other estimations of Ut (see ref- 
erences in the Letter) . Therefore, results for U = At only 
are not sufficient to determine the position of the BOW 
phase with respect to the tricritical point, which is the 
most important qualitative difference between the phase 
diagram in the Letter and those described in Refs. 
To prove the existence of a BOW phase at couplings 
weaker than the tricritical point, one should use param- 
eters U smaller than any estimation of Ut- 

Second, the finite-size-scaling analysis of the charge 
susceptibility Xc(l) m Fig- 1(a) of the Comment is mis- 
leading. A correct analysis is to take the limit N — > oo 
first and then look at the q — > limit. Sandvik, Sen- 
gupta, and Campbell takes both limits simultaneously 
(q = 2tt/N), which can lead to incorrect results. For in- 
stance, the function Fjy(g) = l/(qN) vanishes if the limit 
N — > oo is taken first, but tends to a constant 1/2-k if 
both limits arc taken simultaneously. Thus, the results 
shown in the Comment are no proof of a continuous phase 
transition as a function of V for U = At. 

Third, although I can not rigorously exclude the ex- 
istence of an extended BOW region in the phase dia- 
gram, my results show that its width would certainly be 
much smaller than predicted in Ref. |2j. The main fea- 
tures of the BOW phase (as compared to the competing 
Mott insulator phase) are (i) a long-range-ordered BOW 
(dimcrization) and (ii) a spin gap. I have found a vanish- 
ing spin gap in the thermodynamic limit for the example 
presented in Fig. 1(b) of the Comment. In their previ- 
ous work 0, Sengupta, Sandvik, and Campbell did not 
present any conclusive evidence for the opening of a spin 
gap in an extended region outside the charge-density- 
wave (CDW) regime. It is possible that the spin gap is 
too small to be detected in the finite systems investigated 


(N < 1024 sites), but it is as likely that finite-size ef- 
fects and an arbitrary extrapolation to the infinite system 
limit are responsible for the rather small dimerization re- 
ported in the Comment. I consider that the existence 
of the BOW phase is demonstrated only in those cases 
for which numerical results are consistent. In particu- 
lar, both the extrapolated spin gap and the extrapolated 
dimerization should be clearly larger than zero. 

Fourth, the discrepancies between Sandvik, Sengupta, 
and Campbell results and my results are certainly not a 
failure of the DMRG method nor an effect of open bound- 
ary conditions. In the ground state, the staggered bond 
order of an open finite chain is always larger than in a 
corresponding periodic system because of the Friedel os- 
cillations induced by the chain edges. For both types 
of boundary conditions the staggered bond order ob- 
tained with DMRG decreases with increasing numerical 
accuracy (i.e., an increasing number m of density-matrix 
eigenstates kept). Thus, DMRG results for an open fi- 
nite system systematically overestimate the dimerization 
of the infinite system. The likely cause of the discrepan- 
cies is the difficulty in extrapolating numerical results to 
the thermodynamic limit in the critical region U ~ 2V. 

Finally, the most significant finding in my Letter is the 
presence of the BOW phase at couplings clearly stronger 
than the tricritical point. This fundamentally contradicts 
the theory [2j predicting an extended BOW phase only 
at couplings weaker than (Ut, Vt). Nevertheless, Sandvik, 
Sengupta, and Campbell do not dispute this finding nor 
provide any explanation for this failure of the theory that 
they claim to confirm in their Comment. 

In conclusion, none of the numerical results presented 
in the Comment refute the conclusions of my Letter. 
While the phase diagram presented in the Letter is par- 
tially based on some hypotheses, it is supported by reli- 
able numerical results and a consistent theory. 
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